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Intensive light pulse interaction with a dense resonant medium is considered. The possibihties of 
optical switching and pulse compression at realistic parameters of the medium are analyzed. Pulse 
shape transformation in different photonic band gap structures containing a dense resonant medium 
is studied. In particular the effect of dispersion compensation due to nonlinear interaction with a 
medium is reported. The possibility of pulse control with another pulse is considered in the schemes 
of CO- and counter-propagating pulses. It is shown that a photonic crystal makes controlling more 
effective, at least in the case of co-propagating pulses. 



PACS numbers: 42.65.Re, 42.65.Pc, 42.70.Qs 



I. INTRODUCTION 

The model of resonant two-level medium is studied ac- 
tively for many decades (see, for example, [JQ)- This 
interest is caused, in many respects, by the possibility of 
ultrashort pulses generation which can be used in differ- 
ent fundamental and practical applications. The concept 
of the so-called dense resonant media implies the neces- 
sity of taking into account the interaction of an atom with 
local field produced by all other particles of the system. 
In other words, one has to consider near dipole-dipole 
(NDD) interactions which are characterized by the quan- 
tity 



b = 



(1) 



where /i is the transition dipole moment, C is the volume 
density of two- level atoms, 72 — I/T2 is the rate of trans- 
verse relaxation, h is the Planck constant. Nonlinear ef- 
fects connected with NDD interactions become noticeable 
when the medium is dense enough. One of them, the ef- 
fect of intrinsic optical bistability (JOB), was studied in 
details [1, |3i S @| ■ It was shown that lOB occurs when b 
is greater than 4 (in the case of thin films 0) or even less 
(if we consider propagation effects in extended medium 
Q). This phenomenon takes place in the steady state, 
while in pulse regime one can observe optical switching 
and soliton formation [ll|, [l^] ■ In the present pa- 
per we also consider intensive light pulses, so that change 
of the population difference (or inversion) cannot be ne- 
glected. To study pulse propagation, the numerical so- 
lution of the Maxwell wave equation was implemented. 
This allows to automatically take into consideration the 
processes of dispersion and diffraction which are usually 
ignored in the analytical approaches. 

The main attention in this research is devoted to 
discussion of the properties of a combination of a 



one-dimensional photonic crystal and a dense resonant 
medium under pulse operation. The stationary bistable 
response of such system was analyzed previously [8] . It 
turned out that one can change and control such charac- 
teristics of bistability as hysteresis loop width and switch- 
ing intensity by using photonic crystal. In the present 
paper pulse form transformation in such nonlinear pho- 
tonic band gap structures is studied in comparison with 
the behavior in linear case. In particular, possibilities to 
compensate dispersion spreading in such system and con- 
trol pulse properties by using another pulse are expected 
to be found. 

The paper is divided in several sections. In Section HIl 
the main expressions for description of the system con- 
sidered are given, just as brief characterization of the 
numerical approach used. This methodology is then ap- 
plied to obtain the results of the other sections. Section 
mil is devoted to optical switching and pulse shape change 
(compression and splitting) under propagation in a finite 
layer of a dense resonant medium. In particular, the role 
of near dipole-dipole interactions is estimated. Section 
IIVI contains the calculation results for a single pulse in a 
nonlinear photonic crystal, such as possibility of disper- 
sive spreading compensation. Finally, in Section |V] two 
schemes of controlling pulse intensity with another pulse 
are considered. Photonic crystal is studied as an element 
which makes this process more effective, at least in the 
scheme of co-propagating pulses. 
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II. BASIC EQUATIONS AND NUMERICAL 
APPROACH 



Let us consider radiation propagation in a dense reso- 
nant medium in z-direction. Light interaction with reso- 
nant medium with taking into account nonlinear contri- 
bution due to near dipole-dipole interactions is described 
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by the modified Maxwell-Bloch system as follows [13lll4 



dP 
Itt 
dN 
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where N is the population difference, P is the microscopic 
(atomic) polarization; Auj is the detuning of the field 
frequency to from atomic resonance; 71 = 1/Ti is the rate 
of longitudinal relaxation; c is the light speed in vacuum. 
Macroscopic electric field S is expressed via its amplitude 
iJ as E = i?exp[— i(a;i — kz)]; similarly for macroscopic 
nonlinear polarization we have Pni = fJ-C P exp[—i{ujt — 
kz)]. Here fc = w/c is the wavenumber, and etg is the 
background dielectric permittivity, assumed to be linear 
and dispersionless. 

The system ([^M]) can be represented in the dimen- 
sionless form by introducing new arguments t = cot and 
^ = kz: 



dP 
dN 
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where fl = {fi/hLu)E is the dimensionless amplitude of 
electric field (normalized Rabi frequency); S = Alo/lu 
is the normalized frequency detuning; e = Anfi^C /3huj — 
72& is the NDD interaction constant which provides extra 
nonlinearity in Eq. ((5|); = 7j/w, j = 1, 2. 

Numerical solving of the system ([5]l7]) is performed us- 
ing the finite-difference time-domain (FDTD) method. 
The explicit scheme to solve Eq. ([7]) on the computa- 
tional mesh {I At, jA^) is given by 



where 



ai = ej{l + iAr), 
2 
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Ej is the value of dielectric permittivity Sbg at = j A^. 
In Eq. ([5]) R stands for the finite-difference representa- 
tion of the right-hand member of Eq. ([7]) which describes 
nonlinear properties of the medium. For the case of the 
dense resonant medium it is 

R^3e {P'+\l - iAr) + Pj-i(l + iAr) - Pj(2 + Ar^)) . 




o 

CL 



0.000 0.005 0.010 0.015 0.020 

Normalized amplitude 

FIG. 1: (Color online) Population difference on the entrance 
of the layer of the dense resonant medium after pulse passage 
versus its amplitude Qq at different values of NDD interaction 
constant. Layer thickness L — X. 



Polarization in the mesh points is obtained from Eqs. 
([5]) and ([6]). To solve them the well-known midpoint 
trapezoidal method is used. To set the boundary con- 
ditions we apply the TF/SF approach, when the full cal- 
culation region is divided into two subregions containing 
total and scattered fields respectively, and the so-called 
absorbing boundary conditions according to the perfectly 
matched layer (PML) method. This approach helps to 
avoid nonphysical refiections of the scattered field back 
into the calculation region p^ . 

The calculation scheme considered is similar to that 
of Ref. [l^ where it was used to analyze validity of 
the quasiadiabatic approximation for consideration of 
ultrashort light pulses propagation in a dense resonant 
medium. 



III. COHERENT PULSE INTERACTION WITH 
A DENSE RESONANT MEDIUM 

As it is known, the cases of coherent and incoherent 
interactions are distinguished by the comparison of the 
pulse duration tp with the characteristic relaxation times 
of the medium Ti and T2, moreover, as a rule, Ti > T2. 
Here we consider the coherent case when tp << T2 < Ti . 
In this section the values of parameters are assumed to 
be as follows: light wavelength A = 0.5 /im, tp — 30 fs, 
Ti = 1000 ps, T2 = 100 ps, 6 = 0, and, if not stated an- 
other, Ebg = 1. The pulse on the entrance of the medium 
has Gaussian shape E = Eo exp{—t' /2tp). Obviously, in 
time intervals comparable with the pulse duration (of the 
order of 0.1 ps) the processes of incoherent relaxation can 
be neglected. Under these conditions the system consist- 
ing of two-level atoms can demonstrate optical switching 
from the ground state to the excited one for the time of 
the order of the pulse duration 0. However, as the au- 
thors of Ref. [l3| note, taking into account for the prop- 
agation effects (hence, self-phase modulation as well) not 
only changes the characteristics of switching, but also 
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FIG. 2: (Color online) Calculation results for pulse propaga- 
tion through the layer of the dense resonant medium. The 
amplitude of pulse is Qq = I.SAQt- The thickness of the 
layer L — X; NDD interaction parameter b = 2000. 



makes it difficult to qualitatively predict them without 
execution of rigorous numerical simulations in every par- 
ticular case. 

The calculation results in Fig. [1] show that, at small 
values of the NDD interaction constant 6, the dependance 
of population difference on the pulse amplitude has pe- 
riodic form. When the parameter of extra nonlinearity b 
is increased, this strictly periodic situation is disturbed, 
inversion maximum being shifted towards greater ampli- 
tudes and even being reduced, requiring more accurate 
adjustment of pulse intensity. In addition, there is sim- 
ilar dependance on the layer thickness: strictly periodic 
behavior of population difference occurs only at small (as 
compared with wavelength) medium thickness. 

The difficulties of optical switching, noted in Ref. p^ . 
are connected with the features considered above, as they 
have analyzed the case of large values of NDD interaction 
constant. In our designation it corresponds to the values 
of b of the order of hundreds and thousands. It seems not 
to be realistic: usually b does not exceed several units. 
For example, for gaseous media with typical parameters 

^2 ^ iQ-38 gj.g ^^3^ ^2 ^ iq9 g-l^ (J ^ j^o^O cm-3 [|, [3 

one obtains 6 « 4. In the case of excitonic media pos- 
sessing substantially greater dipole moment (u^ = 10~^^ 
erg cm^, 72 = 10^*' s^^, C ^ 10^^ cm'^ [l||) we have 
only b « 0.4. The dramatic increase of the parameter of 
NDD interaction is unlikely to expect. In the present pa- 
per we keep to this restriction, so the periodic behavior 
corresponding to small values of b will be always valid. 

This implies that the effect of NDD interaction on 
pulse propagation in a dense resonant medium is not sig- 
nificant at realistic parameters. To directly ascertain this 
conclusion, one should compare the results of the calcu- 
lations, when the term including parameter b is present 
or absent in Eq. ([2]) or ([5]). It turned out that the results 
remain the same in both cases. It seems to be obvious, if 
we recall that the parameter b amounts to several units 
(at best, b = 10). This corresponds to the utterly small 
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FIG. 3: (Color online) Temporal behavior of population dif- 
ference on the entrance of the layer of the dense resonant 
medium at the amplitude of the pulse (a) Qq — nAilr, (b) 
Qo — {n ~ 1/2)AQ,T- Parameter 6=1. 



1.0-1 




FIG. 4: (Color online) Temporal behavior of population dif- 
ference on the entrance of the layer of the dense resonant 
medium at different values of the parameter b. The ampli- 
tude of the pulse is f2o = 1.5AS7t. 



value of e w 10^^ — 10"^ in normalized equation ([5]). Only 
when the NDD interaction constant takes on sufficiently 
large values, the contribution of this term becomes es- 
sential (Fig. [21). Thus, while in stationary regime NDD 
interaction plays vital role in such nonlinear phenomena 
as intrinsic optical bistability Q , for coherent pulse prop- 
agation this influence seems to be negligible, at least for 
realistic values of parameters. 
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FIG. 5: (Color online) Distribution of population difference 
along the thickness L = A of the layer of the dense resonant 
medium at different values of the NDD interaction parameter. 
The amplitude of the pulse is flo = LSA^t- 



Since the influence of NDD interactions is negligible, 
one can estimate the period AVIt of the dependance 
shown in Fig._^ by using the commonly known concept 
of pulse area 



Edt. 



(9) 



The period Afiy corresponds to the pulse which returns 
atoms after excitation exactly into the ground state. The 
area of such pulse should be equal to 2tt. Thus, for Gaus- 
sian pulse we obtain 



A 



2v 27rcii] 



(10) 



For the case considered (Fig. [T]) it is approximately 
AUt ~ 0.011. If amplitude of incident pulse amounts 
to an integer number of these periods, i.e. = uAVLt, 
n = 1,2,3,..., then in temporal behavior of population 
difference (Fig. [3^) one can observe n minima and n — 1 
maxima before it achieves a stationary level, in this in- 
stance corresponding to the ground state of the system 
(iV = 1). In turn, at flo = {n — l/2)Af2T n — 1 maxima 
and minima occur (Fig. [8)3) before the medium becomes 
inverted {N = — 1). However, this final value of popu- 
lation difference after pulse passage is not stable, as far 
as relaxation to the ground state becomes apparent at 
longer time intervals. It proceeds the faster, the greater 
constant b (Fig. [5]). Obviously, incoherent relaxation 
does not have enough time to appear. Therefore one 
can attribute this effect to nonlinear interaction of radi- 
ation with the resonant medium. Moreover, this relax- 
ation process should be apparent in spatial scale giving 
rise to nonuniform distribution of population difference 
along the layer thickness (Fig. The similar behavior 
also takes place at smaller values of b, but significantly 
larger thicknesses of the layer. 

Since the infiuence of NDD interactions is negligible, 
the formation of stationary pulses (solitons) of sech-shape 
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FIG. 6: Pulse form transformation after propagation through 
the layer of the dense resonant medium of different thick- 
nesses. The amplitude of pulse is Qq = l.bAflr- NDD inter- 
action parameter b — 10. 
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FIG. 7: Pulse form transformation after propagation through 
the layer of the dense resonant medium of different thick- 
nesses. The amplitude of pulse is fio = 2Af2T. NDD interac- 
tion parameter b — 10. 
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FIG. 8: (Color online) Reflection spectra of photonic crystal 
(number of periods 8) for different values of the thickness d2: 
(a) without a defect, (b) with a defect of L = 20A and — 1. 



can be expected. At the same time energy conservation, 
area change and shape transformation (from Gaussian to 
hyperbolic secant) may cause the pulse to become shorter 
and more intensive. However, this process is limited by 
diffraction and dispersion of light: pulse begins to spread 
as it propagates in the medium. This fact is the reason of 
occurrence of the length of optimal pulse compression. In 
Fig. [S]it amounts about 120A. Note that, after reaching 
this first maximum, several other are observed: at L w 
450A (Fig. [6b), i w 700A and L « 900A (Fig. ifc). 
Finally, at certain (sufficiently large) distance the pulse 
exhibits ultimate attenuation and decay (Fig. [B]d). 

Another important effect that can be observed in this 
system is pulse splitting into several components. The 
number of components depends on the area of initial 
pulse. For example, if the amplitude of the pulse is 
rio = SAi^T (area equals Att), it undergoes splitting into 
two ones (Fig. [7]) corresponding to two 27r-pulses. The 
second (low-intensive) pulse retards more and more from 
the first one as they propagate in the medium. Simulta- 
neously both these pulses lose energy with distance due 
to diffraction and dispersion. 



IV. SINGLE PULSE IN PHOTONIC CRYSTAL 
WITH A DENSE RESONANT MEDIUM 

As it was already stated, the changes of laser pulse 
characteristics under propagation in nonlinear photonic 
band gap (PBG) structures are of special interest in the 
present paper. The one-dimensional photonic crystal to 
be considered here is a sequence of two alternating lay- 
ers with different thicknesses and background refractive 
indices (further they are assumed to equal ni = 1 and 
712 = 3.5). The thickness of the first layer is di = 0.4 /im, 
while the second one, d2, is treated as a variable parame- 
ter which allows to change spectral characteristics of the 
periodic structure. The reflection spectra for different 




without defect 

nonlinear defect 
linear defect 




FIG. 9: (Color online) The forms of transmitted (a, c, e) and 
reflected (b, d, f) pulses after interaction with the photonic 
crystal (number of periods 8) with a defect layer of thickness 
L = 20A, n3 = 1, 6 = 10. The thickness d2 is variable: (a, b) 
da = 0.125, (c, d) da = 0.127, (e, f) da = 0.13 fim. 



values of d2 are shown in Fig. [8^. They are plotted in 
the vicinity of the main wavelength A = 0.5 /im. The 
pulse amplitude in this section is considered to be equal 

to 1^0 = l-^AflT- 

Light pulse propagation in photonic crystals, possess- 
ing strong dispersion due to periodic variation of optical 
properties, results in their spreading on large distances. 
As it is known [13, [3, for materials with sufficiently 
strong nonlinearity, the competition between dispersion 
and nonlinear interaction leads to effective compression 
of pulses. These results regard media with Kerr-type 
nonlinearities. Here we consider some properties of the 
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FIG. 10: (Color online) The forms of transmitted pulses after 
interaction with the photonic crystal with linear and nonlinear 
layers d\ (n\ = 1, 6 = 10). The number of periods is (a) 100, 
(b, c, d) 250. The thickness is variable: (a, b) d2 = 0.13, 
(c) di = 0.127, (d) da = 0.1285 /.im. 
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FIG. 11: (Color online) Reflection spectra of photonic crystal 
(number of periods 250) for different values of the thickness 

d2. 
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FIG. 12: (Color online) The forms of transmitted pulses after 
interaction with the photonic crystal (number of periods 250) 
with both nonlinear layers (n\ = 1, na = 3.5, h = 10) for 
different thicknesses d^. 



systems combining PBG structure with the dense reso- 
nant medium (the value of the NDD interaction constant 
is 5= 10). 

The first such system is a structure containing a defect 
layer in the middle of a sequence of alternating layers. 
The reflection spectrum of such system with the linear 
defect (the refractive index 71,3 is unity) is shown in Fig. 
|8}d. The well-known feature of the photonic crystals with 
defects is seen: appearance of the defect modes, i.e. nar- 
row spectral peaks in the band-gap. Fig. [5] demonstrates 
the computation results of pulse transmission and reflec- 
tion for linear and nonlinear defect layers. As transmis- 
sion of the PBG structure increases (as ^2 changes), the 
effect of nonlinear interaction on the pulse shape becomes 
stronger leading to more effective compression (Fig. [3]3). 
As to reflected radiation, there are several pulses, the first 
one being reflected directly from the photonic crystal lay- 
ers. The others were reflected after interaction with the 
dense resonant medium. 

It seems to be more interesting to consider a vari- 
ant when the layers of photonic crystal are nonlinear. 
Fig. [TO] gives the simulation results when the layers d\ 



of the PBG structure are flUed with the dense resonant 
medium. It is seen that use of nonlinear layers allows to 
effectively compensate dispersive spreading of pulse oc- 
curring in the case of linear photonic crystal. Different 
values of the layer thickness di correspond to different 
values of reflectivity of the linear structure for the ra- 
diation wavelength A = 0.5 ^m (Fig. [TT|) . Naturally, 
more intensive transmitted pulse is observed for smaller 
reflectivity (when = 0.13). At the same time, position 
of wavelength with respect to PBG spectrum influences 
on the duration of pulse transmission through the sys- 
tem: the smaller reflectivity of the structure considered, 
the shorter time interval required for the pulse to appear 
on its output. Comparing transmission durations in Fig. 
llOb -c. one can note that radiation needs only about 67ip 
for d2 = 0.13 /im, while for d2 — 0.1285 /im (the largest 
reflectivity) this pulse delay is equal to 85tp. 

Finally, we consider another system when both layers 
of photonic crystal are nonlinear. On the whole, the be- 
havior of the system (Fig. \Vl\ remains approximately 
the same as in previous case. The maximal change oc- 
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FIG. 13: (Color online) Pulse forms (a-d) and corresponding 
population difference dynamics on the entrance of the medium 
(e-h). Intensity of the first pulse (b, f) fli = 0.5, (c, g) f2i = 1, 
(d, h) Qi — 1.5, while for the second one 02 — 1. Pictures (a, 
e) are for the case of single pulse f22 = 1. Other parameters: 
L = lOOA, b = 10. 



FIG. 14: (Color online) Pulse forms (a-d) and correspond- 
ing population difference dynamics on the entrance of the 
medium (e-h). Intensity of the first pulse (b, f) Qi = 0.5, (c, 
g) Qi — 1, (d, h) Qi = 1.5, while for the second one O2 = 1.5. 
Pictures (a, e) are for the case of single pulse Q2 = 1.5. Other 
parameters: L = lOOA, b — 10. 



curs for the pulse at d2 = 0.1285 fim: the peak intensity 
decreases, while the time delay increases. Note that at 
d2 ~ 0.13 //m the effective splitting of the pulse is ob- 
served. It may be connected with the change of reflec- 
tion characteristics of the nonlinear PBG structure under 
consideration. 

Thus, interaction of light pulse with nonlinear PBG 
structure leads to effective compensation of dispersive 
spreading. On the other hand, using of photonic crystals 
allows to control intensity and time retardation of trans- 
mitted pulse. Additional possibilities to control pulse 
properties are connected with use of a second one. 
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CONTROLLING PULSE INTENSITY BY 
USING ANOTHER PULSE 



FIG. 15: Pulse forms at different distances in medium. In- 
tensities of the pulses Oi — 0.5, Q2 = 1.5. Pictures (b, d) are 
for the case of single pulse O2 = 1.5. Parameter b = 10. 



A. Co-propagating pulses 

Let us consider "one- by-one" propagation of two short 
pulses and interaction between them via the dense reso- 
nant medium. As we have seen in Section IIIIl the pop- 
ulation difference depends on the amplitude of Gaussian 
pulse in periodic manner at realistic values of parameter 



b. In this Section amplitudes will be expressed in the 
units of this period Aflx- Hence, the behavior of the 
second pulse with amplitude O2 differs according to am- 
plitude of the first one fii and, consequently, the state of 
the medium after it has passed. The results of calcula- 
tions are shown in Figs. [01 and [Til Obviously, if the first 
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FIG. 16: (Color online) The forms of transmitted pulses af- 
ter interaction with the photonic crystal (number of periods 
250) with (a) nonlinear layers di, (b) both nonUnear layers. 
Intensity of the second pulse is Q2 ~ 1.5; fli is variable. Pa- 
rameters: ni = 1, n2 = 3.5, di — 0.4 /im, d2 = 0.13 /im, 
b = 10. 
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FIG. 17: (Color online) The forms of transmitted pulse i}2 
controlled with a counter-propagating pulse Sli. Intensity of 
the second pulse is (a, c) Q2 = 1, (b, d) fl2 — 1.5; fli is 
variable. Other parameters: L — lOOA, 6 — 10. 



pulse returns the medium to the ground state (ili = 1), 
it practically does not effect on the second one. But if 
after the passage of the first pulse the medium is excited, 
the second pulse can demonstrate significant increasing 
(Fig. [T3b) or decreasing (Fig. [Mb . d) of peak intensity. 
One can treat this effect as controlling of the second pulse 
intensity by using the first pulse. 

The thickness of the dense resonant medium turned 
out to be an important parameter (see Fig. [TS]) . When 
the distance traveled by the pulse ^2 = 1.5 is L = 50 A or 
L = lOOA, the peak intensity of it appears to be greater 
in the case of the single pulse than in two-pulse scheme. 
But at L = 150 A (when the single pulse demonstrates 
attenuation due to diffraction and dispersion) situation 
becomes reverse. In general, it seems that the first pulse 




FIG. 18: (Color online) The forms of transmitted pulse ^2 
after interaction with the photonic crystal (number of periods 
250) with (a) nonlinear layers di, (b) both nonlinear layers. 
Intensity fl2 = 1.5; intensity of the counter-propagating pulse 
Qi is variable. Parameters: ni = 1, n2 = 3.5, di — 0.4 /im, 
d2 = 0.13 fim, b = 10. 



changes the distance of optimal compression of the sec- 
ond one. 

In order to make the efficiency of control higher, the 
PBG structure can be used. For nonlinear layers di , more 
effective compression is obtained only when fli = 0.5 
(Fig. [TBk). However, when both layers di and d2 are 
nonlinear. Fig. 116b demonstrates increasing of peak in- 
tensities for all values of Sli. For example, for fli — 1.5 
it reaches approximately twofold growth in comparison 
with the single pulse case. The reason of this effect is 
that photonic crystal provides intense energy exchange 
between the first and second pulses due to reflections on 
the layers boundaries. The side effect is the difficulty of 
pulse separation on the output of the system. 



B. Counter-propagating pulses 

Another scheme of controlling pulse intensity is con- 
nected with utilizing of counter-propagating pulse. The 
advantage of this scheme is the convenience of separation 
of pulses as they propagate in opposite directions. The 
calculation results for the pulse ^l2, which is controlled 
with the counter- propagating one fii, are shown in Fig. 
[TTI It is seen that change of the intensity of the incident 
pulse rii leads to the change of the intensity of the trans- 
mitted pulse fl2 , in particular to its decreasing in compar- 
ison with single pulse case. Since the change of the length 
of optimal compression occurs in this scheme as well, one 
can obtain increasing of the intensity of ^2 = 1-5 pulse 
on the distance of about 200 A when 17i = 1. 

Fig. [18] shows that the interaction between pulses in 
the nonlinear photonic crystal is negligible, in contrast 
to the case of co-propagating pulses. It seems to be the 
result of weak coupling of the pulses in such system, as 
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they propagate almost independently. This is the great 
disadvantage in the view of controlling possibilities, but 
it can be used for simultaneous work with two pulses 
moving in opposite directions. 

VI. CONCLUSION 

The consideration of intensive radiation interaction 
with the dense resonant medium in the case of coher- 
ent pulse regime allows to conclude that the influence of 
near dipolc-dipole interactions can be neglected at real- 
istic parameters of the medium, in contrast to stationary 
case. At the same time, one can observe pulse compres- 
sion (and splitting) which has certain optimal distance 



due to the processes of diffraction and dispersion. This 
property of nonlinear compression can be used to com- 
pensate dispersive spreading in photonic band gap struc- 
ture containing the resonant medium considered. 

It turned out that dense resonant medium allows to 
control pulse with another pulse in schemes of co- and 
counter-propagating pulses. It seems to be useful from 
the point of view of prospective techniques of opti- 
cal information storage and processing. Photonic crys- 
tal makes the process of controlling more effective, at 
least in the case of co-propagating pulses. For counter- 
propagating ones, it appears not to be efScient, neverthe- 
less independence of pulse propagation in opposite direc- 
tions can be used for parallel work with two pulses. 
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